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The integrable structure ol Ginibre's orthogonal ensemble of random matrices is looked at through 
the prism of the probability p n ,k to find exactly k real eigenvalues in the spectrum of an n x n 
real asymmetric Gaussian random matrix. The exact solution for the probability function p n ^ 
is presented, and its remarkable connection to the theory of symmetric functions is revealed. An 
extension of the Dyson integration theorem is a key ingredient of the theory presented. 
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Introduction. — In the mid-1960's, Ginibre introduced 
statistical ensembles of (i) real, (ii) complex, and 
(iii) real quaternion random matrices whose eigenvalues 
may belong to any point of the complex plane C. They 
were derived from the celebrated Gaussian orthogonal 
(GOE), Gaussian unitary (GUE), and Gaussian symplec- 
tic (GSE) random matrix ensembles Q in a purely formal 
way by dropping the Hcrmiticity constraint and, there- 
fore, can be thought of as their direct non-Hermitean 
descendants. Respectively coined as GinOE, GinlJE, 
and GinSE, non-Hermitean random matrix models ex- 
hibited intriguingly rich mathematical structures whose 
complexity exceeded by far that of Hermitean random 
matrix theory (RMT). 

From the physical point of view, non-Hermitean ran- 
dom matrices have proven to be as important Q as their 
Hermitean counterparts [4|. Ginibre's random matrices 
appear in the description of dissipative quantum maps |f| 
(GinUE), in the studies of dynamics and the synchro- 
nisation effect in random networks (GinOE), in the 
statistical analysis of cross-hemisphere correlation matrix 
of the cortical electric activity |8| (GinOE), and in the 
characterisation of two-dimensional random space-filling 
cellular structures [9j (GinUE). They also arise in the 
context of "directed quantum chaos" HUEI (GinOE, 
GinUE, GinSE). Chiral deformations of non-Hermitean 
random matrices (GinOE, GinUE, GinSE) help elucidate 
universal aspects of the phenomenon of spontaneous chi- 
ral symmetry breaking in quantum chromodynamics |l2| : 
the presence or absence of real eigenvalues singles out dif- 
ferent breaking patterns. More recent findings 0| asso- 
ciate statistical models of non-Hermitean normal random 
matrices with integrable structures of conformal maps 
and interface dynamics at both classical 14] and quan- 
tum scales [l f5l | . For a comprehensive review of these and 
other physical applications, the reader is referred to Ref. 

n 

Out of the three non- Hermitean matrix models 1], 
Ginibre's orthogonal ensemble defined by the probabil- 
ity density P G i„o E [H] - (2tt)-" 2 / 2 exp [- tr (HH T /2)] 
for an n x ii real matrix TL to occur is the least stud- 
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ied and the most challenging. Perhaps, the great dif- 
ficulties faced in statistical analysis of the GinOE can 
be attributed to the fact that its, generically complex, 
spectrum (w\, ■ ■ • ,w n ) may contain a finite fraction of 
real eigenvalues. This very peculiar feature of GinOE 
can conveniently be accommodated by dividing the entire 
space T(n) spanned by all real n x n matrices Tt £ T(n) 
into (n + 1) mutually exclusive sectors T(n/k) asso- 
ciated with the matrices having exactly k real eigen- 
values, such that T(n) — Ul-=o ^( n /^)- The sectors 
T(n/k), characterised by the partial probability densi- 
ties P-Hei(n/k){ w u * T w n), can be explored separately 
because they contribute additively to the joint probabil- 
ity density function ( JPDF) of all n eigenvalues of TC from 
T(n): 



P n (wi, ■ ■ ■ ,w„) 



y] p-h€T(i 

k=0 



,/k)(wi, ■ ■ ■ ,W n ). (1) 



Two spectral characteristics are of particular physical in- 
terest: (i) the partial (r±, r2)-point correlation functions 
obtained by integrating out all but r\ real and T2 complex 
eigenvalues in Pu&T(n/k) > and (ii) the probability 
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to find exactly k real eigenvalues in the spectrum of 
GinOE that additionally contains £ pairs of complex con- 
jugated eigenvalues, so that n = k + 21. The probability 
function p n jz is the central object of our study, which, if 
considered in a wider context, aims to highlight an exclu- 
sive role played by symmetric functions in the description 
of GinOE. 

The most important and general result available for 
GinOE is due to the breakthrough work by Lehmann and 
Sommers [l6| who proved, a quarter of a century after 
Ginibre's work, that the fe-th partial JPDF (0 < k < n) 
equals [HI 
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In writing © , we have used a representation due to Edel- 
man who rediscovered the result 0] a few years 
later. The above JPDF is supported for (Ai, • • • , Afc) S 
R k , (Rezi, • • • ,Rez e ) e and (Im<?i, • • • ,lmze) <E 
(M+)^. A particular case k — n of ©, corresponding 
to the matrices H £ T(n/n) with all eigenvalues real, 
was first derived by Ginibre 

Equations and © solve the problem of finding the 
JPDF of all n eigenvalues in GinOE. Despite this tremen- 
dous progress, a calculation of the (ri, r2)-point correla- 
tion function based on turns out to be the problem, 
particularly, because the well-developed machinery of 
RMT, which was at use in studies of complex eigenvalue 
correlations in GinUE [l| and GinSE |ig, |20(, fails to 
work in this case. In particular, the famous Dyson in- 
tegration theorem |2l| (Theorem 5.1.4 in Ref. is no 
longer applicable to the study of spectral statistics in 
GinOE, as will be reasoned below. This failure clearly 
signals of novel mathematical structures lurking behind 

To unveil and explore these structures, we wish to con- 
centrate on the probability function Pn.k- Previous at- 
tempts by Edelman and co-workers |18l |22j to attack the 
problem brought no explicit formula for p n ,k for generic 
k. The analytic results currently available include: (i) 
the probability of having all n eigenvalues real equals 
[3 Pn n = 2 -»(«- 1 )/ 4 (this is the smallest probability 
out of all p n ,k)', (h) for all < k < n, the p n ^ is of the 
form p nyk = r n , k + s, hk y/2, where r n> k and s Ui k are 
rational; (iii) the expected number E n = Y^k=o ^Pn,k of 
real eigenvalues is given by j2^| 



-1/2; n; 1/2) 



B(n,l/2) 



(4) 



Main results. — The following four statements sum- 
marise the new results reported in this Letter. 

Statement 1. The probability p n> u of exactly k real 
eigenvalues occurring equals 



Pn,k = Pn,n-1i = Pn,n Fe(Pl, ' ' ' ,Pt)- 

The universal multivariate polynomials 



,Pe) 



M=t 3 = 1 



(5) 



(6) 



are exclusively determined by the partitions A of 
the size |A| = £, where I is the number of pairs 
of complex-conjugated eigenvalues. The formula © 
uses a frequency representation |23| of those partitions 
A = (II 1 , ■ ■ ■ ,£g 9 )- Experts in the theory of symmetric 
functions may readily recognise that our T^s are, up to 
a factorial, so-called zonal polynomials |24j 
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fe(pi, ■■■ ,Pt) = jj Z {lt) {p 



> Pit- 



They are tabulated in the manuscript |25| by Jack and 
can also be efficiently calculated [24| by recursion. On the 
contrary, the arguments pj 's of the multivariate polyno- 
mials in (JjjJ are matrix model dependent: 

Pj = tr( ,[„/2]-i) Q 3 ■ (8) 

(The notation [x] stands for the integer part of x). The 
non-universal matrix g is sensitive to the parity of n. For 
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Exact Solution 



Appr. Numerics 



29930323227453-20772686238032^2 
17592186044416 



0.031452 0.031683 



3(1899624551312^-2060941421503) 
4398046511104 



0.426689 0.427670 



3(2079282320189-505722262348^2) 



8796093022208 



0.465235 0.464098 



252911550974^-27511352125 
4398046511104 



0.075070 0.075021 



15(1834091507-10083960y / 2) 
17592186044416 



0.001552 0.001526 



3(1260495y / 2-512) 
2199023255552 



0.000002 0.000002 



8589934592 



0.000000 0.000000 



TABLE I: Exact solution for p\2,k (second and third column) 
compared to numeric simulations (fourth column) performed 
by direct diagonalisation of 1, 000, 000 of 12 x 12 matrices. 
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n = 2m even, its entries are 

QZ7 = / dyy^-^ev 2 eifc(yV2) (2a + 1) 
Jo L 

x Lf£^-\- 2y 2 ) + 2y 2 L^T» +1 (-2 y 2 )l ,(9) 



while for n = 2m + 1 odd, 



"Odd 

0a, /3 



eZJ - (-4) 



ra—j3 



(2m)! 13 



(2/3)! 

-oven. 



(10) 



Here, L^{x) are generalised Laguerre polynomials. 

We wish to stress that zonal polynomials are the 
integral part of the final solution © , and not merely the 
means of its derivation. This highlights a strong con- 
nection between integrable structure of GinOE, integer 
partitions JfjJ, and the theory of symmetric functions J7J. 

Statement 2. The entire generating function for the 
probabilities p n , n -2i can be reconstructed from ©. The 
fairly compact result 

[n/2] 

Gn{z) = ^ zt Vn,n-2i = Pn,n det(l[„/ 2 ] + Z g) (11) 
1=0 

with the g of needed parity provides us with yet another 
way of computing the entire set of p n , k 's at once. For 
comparison of our analytic predictions with numeric sim- 
ulations, see the Table I. 

Statement 3. The previously unknown JPDF of I 
pairs of complex conjugated eigenvalues of a matrix 
H G T(n/k) is given by the formula (fH)|) . 

Statement 4- This statement is formulated in the form 
of the Pfaffian integration theorem (jl7|l below. (It is 
highly likely that this theorem will have implications far 
beyond the scope of the present Letter). 

Sketch of the derivation. — Because of the space limita- 
tions, only key points of the derivation will be presented 



below; the technical details will be reported elsewhere 

To determine the probability function p n , k , defined by 
(J2J), we first integrate out all (Ai,--- ,\k) in (@ which 
can be viewed as a hybrid of GOE [2[ and GinSE |19ll20l| . 
The most efficient way to do so is to spot that the 
first line in (J3J), if taken without any prefactors therein, 
is proportional to the average characteristic polynomial 
(11.7=1 det(«j — H) det(£j — H)) of a k x k random ma- 
trix H drawn from the GOE with the weight function 
exp(— tr H 2 /2). It admits the Pfaffian representation 



fc! 



nl A 2e ({z,z}) 



pf 



( Z{ , Zj ) K n ( %i 7 Zj ) 
K-n ( %i •> %j ) ( %i i ) 



(12) 



2£x2£ 



with the constant c k = 2 fc / 2 fc! rij=i r(j/2) given by a 
GOE Selberg's integral Q- The Vandermonde deter- 
minant A2f({z, z}) is calculated on the set of complex 
eigenvalues (z±,- ■ ■ ,Z£,z~i,--- ,zg). The kernel function 
K n (z,z') in III is a so-called L»-part El of the GOE 
matrix kernel |2| that can be expressed |30| in terms of 
arbitrary monic polynomials qj : 



_^ n— 1 

K n (x,y) = q]{x) fi jk q k (y). 



(13) 



j,k=0 



The real antisymmetric matrix {/, is determined by its 
inverse 

= o / dxdye-^ +y2)/2 sgn{y-x)q ] {x)q k {y). 

(14) 

Substituting l|T2|) into 101 completes the A-integration 
therein, bringing a new result for the JPDF of £ pairs 
of complex conjugated eigenvalues of 7i € T(n/fc): 



The structure of Q15I1 mirrors that of the JPDF of all complex eigenvalues in the GinSE (see, e.g., Chap. 15.2 in 
Ref. [2j, and Ref. |20j), possibly triggering one to think that the remaining z-integrations in Q could readily be 
accomplished by virtue of the Dyson integration theorem 0, Sadly, this is not the case because the D-part 

K n (z, z') of the GOE matrix kernel does not obey the projection property in the complex plane 

j da(w) K n (z,w) K n (w,z') ^~K n (z,z') (16) 

with respect to the measure da(w) — erfc[(w — w)/i\/2] exp[—(w 2 + w 2 )/2] d 2 w from (|15(l (compare to the Lemma 
15.2.1 of Ref. 2J). Fortunately, the projection property is not actually necessary to carry out the integrations. We 
were able to prove the following Pfaffian integration theorem |26j . 

Theorem. — Let da(z) be any benign measure on z S C, and the kernel function K n (x,y) be an antisymmetric 



K n \Zi • Zj) K n (zi : Zj) 
■Kn \Zi , Zj ) Kji \Z% , Zj ) 



(15) 



2/x2/ 



4 



function specified by 113fl . Then, the integration formula 



n / da ( z j) p f 

j=l J Zj£C 



fi^n \%i 3 %j ) K n yZi , Zj ) 
-^-n (-^z i %j ) K n i^Zi , Zj ) 



Z(i 



2f,x2£ 



-tr 



(O.n-l) 17 >' 



' 2 



tr 



(O.n-l)' 



(17) 



holds, provided the integrals in the l.h.s. exist. 

The Pfaffian integration formula contains particular 
zonal polynomials |24j whose arguments are traces 

of a 3 with 



n-1 
k=0 



da(z) [q k {z) qp{z) - qp(z) q k (z)} 



(18) 

Operationally, the proof [26( is based on a topological 
interpretation of the permutational expansion of Pfaffian 
in i|17|) combined with term-by-term integration. 

Equipped with this theorem, we are ready to complete 
the z-integrations in (0. The calculation is most eco- 
nomic in the basis where monic polynomials qj 's coincide 
with GOE skew-orthogonal polynomials. For n = 2m 
even, they are given by 



q 2j (x) = 2- 23 H 2j (x), 



2-^- 1 )ji/ 23 _ 1 (x), 
(19) 



so that 



A _ l + (-l) fc 



1 + 
2 1 ~3j\ 



Sj,k- 



(20) 



Substituting l|19|) and l|2(J|) into JTSJ, and taking into ac- 
count (|17|l . we reproduce the announced solution (jSJ) 
@ after a bit lengthy but straightforward calculations. 
The case n = 2m + 1 odd can be treated similarly j2|| 
leading to the same solution but with © replaced by 

It remains to establish the result (|llfl for the generating 
function G n (z). It stems from the summation formula 



T- 

r=Q 



,Pr) 



exp ^(-l) 1 - 1 



p r Z' 



(21) 



well known in the theory of symmetric functions |24j . 
Identifying p r = tr( 0) [ n /2]-i)(5 r , the exponent in (|21|) can 
be evaluated explicitly to confirm (|llfl . 

Conclusion. — To summarise, the exact solution was 
presented for the probability p Ui k to find precisely k real 
eigenvalues in the spectrum of an n x n random matrix 
drawn from GinOE. Expressed in terms of zonal polyno- 
mials, the solution associates the integrable structure of 
GinOE with the theory of symmetric functions. 

Certainly, more work is needed to accomplish the spec- 
tral theory of GinOE. (i) In particular, the large-n be- 
havior 26] of the probability function p n: k should be ex- 
amined when the number k of real eigenvalues does or 



does not scale with the matrix dimension n. The large- 
n formulae of this kind would facilitate a comparison of 
our exact theory with existing experimental and nu- 
merical (see the second paper in Ref. data, (ii) 
The calculation of all (n, r2)-point correlation functions 
for GinOE as defined below is yet another important 
problem to tackle. We believe that progress in this direc- 
tion can be achieved through a proper extension of the 
Pfaffian integration theorem 1171) . 
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